Realizing controllable noise in photonic quantum information channels 
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Controlling the depolarization of light is a long-standing open problem. In recent years, many 
demonstrations have used the polarization of single photons to encode quantum information. The 
depolarization of these photons is equivalent to the decoherence of the quantum information they 
encode. We present schemes for building various depolarizing channels with controlled properties 
using birefringent crystals. Three such schemes are demonstrated and their effects on single photons 
are shown by quantum process tomography to be in good agreement with a theoretical model. 
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PACS numbers: 03.65.Yz, 42.25.Ja, 42.50.Lc 



Light depolarization is a fundamental optical phenom- 
ena. It was studied as early as the nineteenth century, 
when measurement and characterization methods of the 
polarization state of light were introduced [1]. Meth- 
ods for complete depolarization of hght, such as the 
Cornu and Lyot depolarizers, have been known for many 
decades [i-Q. The Cornu and wedge depolarizers require 
the light beam to be wide, as the loss of coherence be- 
tween the two polarizations is achieved via averaging over 
the spatial degrees of freedom. In the case of the Lyot 
depolarizer, short coherence length is required and the 
averaging is over the temporal degrees of freedom. Nev- 
ertheless, there are currently no methods that enable con- 
trol of all of the aspects of the depolarization process. 

Quantum information is commonly encoded in the po- 
larization of single photons[5]. Depolarization of such 
photons acts as quantum noise on the stored informa- 
tion, i.e. the interaction between the information en- 
coding units and the environment results in decoherence. 
In order to study quantum decoherence in general and 
its effect on quantum information protocols in particu- 
lar, it is desirable to create quantum channels with con- 
trolled noise. Such channels will be useful for testing 
quantum error correction and quantum key distribution 
protocols 6, 7.]. Other uses for these channels are to test 
for the existence of decoherence- free subspacesjl] and for 
generating partially- mixed entangled states j^. 

In recent years, depolarizing channels were studied by 
several methods. When a single birefringent crystal was 
used, only dephasing channels were demonstrated, as we 
will show later[8]. Optical scatterers such as emulsions, 
multi-mode fibers and ground glass, give variable depo- 
larization, but are also accompanied by a spread in k 
space, resulting in considerable loss when collected for 
further processing[9l-[llj| . Another approach is to use po- 
larization scramblers of various kinds |12l. [l3| . These real- 
izations are equivalent to fast polarization rotations and 
averaging measurements over times longer than the typi- 
cal rotation periods. Nevertheless, when used with single 
photons, each photon by itself is completely polarized. 
Controllable depolarizers were demonstrated by using 
two wedge depolarizers with variable beam diameter [l^. 



or with a tunable relative angle[15|. These channels are 
hard to model and their anisotropy level is uncontrolled 
as they couple the polarization with many spatial degrees 
of freedom. 

In this letter we present a theoretical and experimen- 
tal study of various controllable depolarizing channels. 
We study channels that are composed of a sequence of 
birefringent crystals and wave plates. The depolariza- 
tion and its anisotropy depends on the order and rela- 
tive angles between the channel components. The gen- 
erated channels are mostly anisotropic and can be tuned 
continuously between no depolarization and complete de- 
phasing. These channels were characterized by the trans- 
mission of polarized single photons, generated by spon- 
taneous parametric down-conversion. Quantum process 
tomography (QPT) was used to compare the experimen- 
tal results with theory (lil, [l7j. 

The information in a classical channel can be degraded 
only by bit-flip errors. Thus, such a channel is completely 
described by a single parameter - the bit-flip error prob- 
ability. In comparison, quantum channels can have a 
constant unitary rotation and three types of errors, rep- 
resented by the Pauli operators: a bit-flip that swaps the 
logical |0) and |1) amplitudes, a phase flip between the 
amplitudes and the combination of the two, which is a 
third orthogonal operation. Isotropic decoherence is the 
case when the three error probabilities are equal. 

A polarization qubit can be described either by a den- 
sity matrix operator p or equivalently, by a point in the 
Poincare sphere. The Cartesian coordinates of this point 
are the Stokes parameters S = {Si, S2, S3) which de- 
scribe the linear {\h),\v)), diagonal {\p) = (|^) + \ v))/V^, 
\m) = {-\h) + \v))/V2) and circular (|r) = {\h) + 
i\v))/V^, \l) = {i\h) + \v)) / ^/2) polarization components, 
respectively. The Degree Of Polarization (DOP) is de- 
fined as D, the length of the Stokes vector [18j: 



D = JS^ + Si + Sl = ^l^Adet{p). 



(1) 



The perfectly polarized states are described by the sur- 
face of the sphere [D = 1) and its center designates the 
completely unpolarized state [D — Q). The inside of the 
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sphere includes all partially polarized states (0 < D < 1). 
The physical meaning of the DOP is the ratio between 
the polarized light intensity and the total light intensity. 

Consider an arbitrarily polarized wave packet that is 
passing through a birefringent crystal[l^. The tempo- 
ral walk-off T = between two wave packets, each 
polarized along one of the symmetry axes of the crystal, 
depends on the crystal length L, its refractive index dif- 
ference Art and the speed of light c. We assume that 
the coherence time of the wave packets tc is shorter than 
the walk-off r. If the light wave packet is not polarized 
linearly along one of the crystal symmetry axes, its two 
components acquire temporal distinguishability. Thus, 
the polarization and temporal degrees of freedom become 
entangled. The role of the environment in general deco- 
herence models is fulfilled here by the temporal degrees 
of freedom. As the detectors are insensitive to the short 
temporal walk-offs, they cannot distinguish between the 
wave packets, effectively tracing out the temporal degrees 
of freedom. The result is an effective depolarization since 
no coherence can be observed between the two orthogo- 
nal polarizations. The depolarization operation is de- 
scribed in the Poincare sphere picture by a projection of 
the initial Stokes vector on the direction that represents 
the crystal principal axes. For example, a birefringent 
crystal aligned along the h-v directions will project any 
initial state onto the Si direction. This kind of operation 
is referred to as a 'dephasing channel' [5|. 

A single crystal configuration can apply any level of 
depolarization to any initial linear polarization. On the 
other hand, for such a configuration there is always an- 
other polarization direction that experiences no depolar- 
ization whatsoever. Thus, we consider a second crystal, 
that is placed after the first one [4]. For historical reasons, 
let us first assume that the second crystal is twice as long 
as the first one. The two crystals are coupling orthogonal 
temporal degrees of freedom, as the first crystal couples 
t = with r, while the second couples t — with 2t and 
t = T with 3t. Thus, this configuration can be described 
as two consecutive projections of the initial polarization 
state onto the Stokes directions defined by the crystals' 
axes. In the case where the orientation of the crystals 
differ by 45°, the two projections are perpendicular, re- 
sulting in a final state at the sphere center (D = 0, the 
completely unpolarized state) for any initial state. This 
configuration in known as the 'Lyot depolarizer' 

The relevant error rates for practical tests of quantum 
information protocols are below 20%[23|- For this reason, 
it is desirable to have a depolarizing scheme that can be 
tuned to small values of depolarization or even to zero 
depolarization. Thus, we investigate a configuration of 
two identical crystals. If the second crystal is oriented at 
90° with respect to the first one, the polarization time 
delay created by the first crystal is exactly compensated 
for by the second. For any other relative angle, there can 
be up to three different temporal modes: t = Q and r are 
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FIG. 1. (Color online) The experimental setup and depolar- 
izing schemes. See full description in text. 



coupled by the first crystal, while the second crystal cou- 
ples between them and an additional third delay t = 2r. 
Changing the relative angle between the two crystals af- 
fects the occupation of each of the three modes, which re- 
sults in a different depolarization. Hence, by tuning this 
angle we control the channel depolarization level. The 
modal occupation can be easily calculated for a given 
initial polarization and tuning angle. The polarization 
two-dimensional density matrix p/ of the final polariza- 
tion state is the normalized sum of the polarization states 
of all the occupied temporal modes lipt)'- 



Pi 



(2) 



where T, t and at are the number of relevant temporal 
modes, their index and weights, respectively. This calcu- 
lation is equivalent to partially tracing over the temporal 
degrees of freedom of the 2T-dimensional polarization- 
time density matrix. In addition, the effect of wave-plates 
is simple to calculate as they rotate each of the temporal 
modes separately. With this method, it is possible to cal- 
culate the final polarization state for every initial state 
that passes through a sequence of birefringent crystals (a 
"depolarizer"). From the resulting density matrix it is 
possible to calculate the DOP by Eq. [TJ 

A 780 nm Ti:Sapphire pulsed laser of 76 MHz repeti- 
tion rate was frequency doubled and the 390 nm pulses 
were focused into and collinearly down-converted in a 
1mm thick type-I BBO crystal (See Fig. [1]). The down- 
converted signal was filtered by a dichroic mirror (DM) 
and collimated with a lens (L2). One photon of the pair 
was split by a beam-splitter (BS) and detected, and the 
second photon was directed to the depolarizer. The po- 
larization state of the depolarized photons was charac- 
terized by wave plates and a polarizer (POL). Photons 
were filtered by 5nm bandpass filters (IF), correspond- 
ing to a coherence time of tc — 180 fs, and then cou- 
pled into single-mode fibers leading to single photon de- 
tectors (DET). We characterized the depolarization of 
three initial states \p), and \r) which are mut ually 
unbiased^, by quantum state tomography (0ST)[l7|. 
The \v) state was also measured as required by QPT. 




FIG. 2. (Color online) Experimentally measured final states 
for \h) (blue squares), |p) (red triangles), and |r) (green cir- 
cles) inputs, plotted for a range of crystal angles up to a final 
value of (a) £» = 45°, (b) 9 = 54.74°, and (c) 6 = 67.5°. The 
wire-mesh ellipsoids represent the mapping of the surface of 
the Poincare sphere by the depolarizing operation as mea- 
sured by QPT. (d) Comparison between measurements in the 
range 0° < < 90° and the theoretical model. Theoretical 
curves are presented as solid lines in the range 0° < 6 < 180° . 



The first depolarizing scheme that we present is com- 
posed of two 2 mm long Calcite crystals (C) with two A/2 
wave plates (See Fig. [T]). The crystals are fixed perpen- 
dicularly, with one wave plate before and the other after 
the first crystal. Rotation of the wave plates in opposite 
directions by an angle of 9/2 is equivalent to the rotation 
of the first crystal by 6. When 6* = 0° no depolarization 
occurs, and when 6 = 90° the depolarizer is equivalent to 
a dephasing channel of a single crystal. Wave plates were 
used instead of direct rotation of the crystal to eliminate 
an unwanted angle dependant retardation. 

The transformation of purely polarized states through 
the depolarizer for various angles between the two crys- 
tals are presented in Fig. [2l In Figs. [2ja-c), the measured 
QPT mappings of the initial sphere are shown for three 
specific cases. These cases are when the polarization of 
only one initial state is completely lost {9 = 45°), when 
the channel is isotropic {9 = 54.7°) and when two initial 
states are depolarized identically {9 = 67.5°). The pro- 
cess fidelities, when compared to the model for all three 
cases, were higher than 97%. The measured output states 
for the three mutually unbiased polarizations were ob- 
tained using the maximal- likelihood QST algorithm 17 1. 
They are also presented, up to the respective angles. Fig- 
ure [2l^d) compares these results with theory. The results 
show a very good agreement with the model. From these 
results, we have calculated the DOP as a function of 9. 
The DOP results are presented in Fig. [3]with their theo- 
retical predictions. The three special cases of Figs, [^fa-c) 
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FIG. 3. (Color online) Experimentally measured degree of po- 
larization of the output states as a function of the equivalent 
first crystal rotation angle 9 for the first depolarizing scheme. 
Initial states and their representation are the same as in Fig. 
(2] Model predictions are presented as solid lines. 



are clearly reproduced in our measurements. 

In general, this depolarizer is anisotropic, except for 
9 = 54.7° = tan~^(\/2). Analytical calculation of the 
DOP for this angle shows that the final value of 1/3 is 
independent of the initial state, as can be seen in Fig. 
|3l where all three curves intersect. The isotropy of this 
configuration is apparent in Fig. EJ^b), where the polar- 
ization sphere is mapped to another sphere of radius 1/3. 

We studied a second depolarizing scheme that was 
composed of two perpendicularly fixed identical crystals 
with a quarter- wave plate between them (See Fig. [T]). 
The quarter- wave plate angle 9 is set to be zero when the 
principal axes of the wave plate and the first depolariz- 
ing crystal are parallel, for any given initial polarization 
(5i,S'2,5'3), the final DOP is 



■ cos (40) 
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^-- + - COS (40) + -COS 



(3) 



All states with the same \Si\ value, result in the same 
DOP for a certain 9. It is possible to find three mu- 
tually unbiased polarization bases that have the same 
S*! = i-^ value, and thus, experience the same depolar- 
ization. We define this situation as symmetric depolar- 
ization. For such bases, the DOP can be tuned between 
1 and ~ 0.41 as a function of 9. We generated these 
states and characterized them after passing through the 
second depolarizing scheme. The DOP results are shown 
in Fig. |4ja), and their Poincare representation in Fig. 
Sljb). A good agreement with theory is observed. 

The third scheme adds the possibility for symmetric 
depolarization down to complete depolarization. The 
difference between the second and third schemes is the 
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FIG. 4. (Color online) (a) Experimentally measured degree of 
polarization as a function of the quarter-wave plate angle 9 for 
the second and the third depolarizing schemes. For the second 
scheme, the three orthogonal initial states with 5*1 = — ^/T/S 
after depolarization are presented by black squares, blue tri- 
angles and red inverted triangles. For the third scheme, the 
depolarized states of all three initial states should be iden- 
tical. Results for one of these states are presented by black 
circles. (b,c) The measured states in the Poincare sphere for 
the respective initial states. 



doubling of the second crystal's thickness (See Fig. [T]). 
As before, the final DOP depends only on the initial Si 
value and the wave plate angle 9, but now it takes values 
between and « 0.58. This result is due to an effec- 
tive additional Si projection to the output of the second 
scheme. At the 9 = 45° position, the third scheme is ex- 
actly a Lyot depolarizer that completely depolarizes any 
initial polarization by consecutive Si and S^ projections. 

Results for an initially polarized state with Si = — \= 

v3 

are shown in Figs. SKa) andlD^c). As can be seen, when 
9 = 45°, the state is completely depolarized. The state 
tomography results (Fig. |4l[c)) reveal the difference from 
the previous scheme as an extra Si projection. 

Although we have demonstrated the effects of the var- 
ious depolarizing schemes using single photons, these de- 
polarizers would also be effective on any classical light 
with short enough coherence time. We repeated our mea- 
surements with laser pulses, and demonstrated identical 
results (not presented here). Thus, these results apply 
not only to polarization encoded qubits, but to any clas- 
sical scenario where controlled depolarization is required. 

In this work we use crystals that are long enough to 
completely separate the two polarization components. It 
is possible to deliberately use shorter crystals that will 
leave a portion of the two wave packets overlapping. The 
same qualitative results will be achieved this way, but 
with smaller magnitude. Another scheme that enables 
control of the channel isotropy level is to combine three 



setups of the same kind. The crystal lengths should be 
tripled for the second setup and tripled again for the 
third, in order for each of them to affect orthogonal tem- 
poral modes. If the second setup is positioned at 45° 
with respect to the first and a quarter-wave plate at 45° 
is placed between the second and the third setups, then 
each setup will be oriented at a different Stokes direc- 
tion. This scheme can generate isotropic channels of any 
depolarization amount. 

In conclusion, we have demonstrated a scheme for the 
realization of various quantum channels for photon po- 
larization qubits with controllable decoherence. Isotropic 
and anisotropic depolarization processes are possible. 
Channels were characterized by QPT using the maximal- 
likelihood algorithm. All the results are in a good agree- 
ment with a simple theoretical model. These depolarizers 
can be used to evaluate the performance of quantum er- 
ror correction and quantum key distribution protocols. 
In addition, they can be utilized in any classical optics 
setup where controllable depolarization is required. The 
authors would like to thank the Israeli Science Founda- 
tion for supporting this work under grant 366/06. 
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